First-principles study of nucleation, growth, and interface structure of Fe/GaAs 
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We use density- functional theory to describe the initial stages of Fe film growth on GaAs(OOl), 
focusing on the interplay between chemistry and magnetism at the interface. Four features appear 
to be generic: (1) At submonolayer coverages, a strong chemical interaction between Fe and sub- 
strate atoms leads to substitutional adsorption and intermixing. (2) For films of several monolayers 
and more, atomically abrupt interfaces are energetically favored. (3) For Fe films over a range 
of thicknesses, both Ga- and As-adlayers dramatically reduce the formation energies of the films, 
suggesting a surfactant-like action. (4) During the first few monolayers of growth, Ga or As atoms 
are likely to be liberated from the interface and diffuse to the Fe film surface. Magnetism plays an 
important auxiliary role for these processes, even in the dilute limit of atomic adsorption. Most of 
the films exhibit ferromagnetic order even at half-monolayer coverage, while certain adlayer-capped 
films show a slight preference for antiferromagnetic order. 
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I. INTRODUCTION 

Magnetic materials are widely anticipated to be inte- 
grated into semicondudjor-based microelectronics during 
the next decade or two.EJ A major component of this effort 
has focused on using ferromagnetic thin fihus as a source 
for creating spin-polarized electrical current in a semi- 
conductor substrate, a process referred to as "spin injec- 
tion." Three criteria have been identified as important 
for technologically useful spin injection: (1) Substantial 
spin polarization of the injected current; (2) Electrons, 
rather than holes, serving as the spin-polarized carriers; 
(3) Curie temperatures for the source of order room tem- 
perature or higher. 

Although several materials appear to meet one or two 
of these criteria, none has yet met all three. Several di- 
luted magnetic SjSiniconductors based on Be-[-.and Mn- 
doping of ZnSeJj'El and Mn-doping of GaAs,cl have re- 
cently demonstrated spin-injection efficiencies of greater 
than 50%, but only near liquid helium temperatures. 
Other magnetic semiconductors, including CdCrSe, can 
be doped both p- and n-iype and have Curie tempera- 
tures of order 100-200 K,lI but their performance as spin 
injection sources has yet to be examined. 

In addition to the magnetic semiconductors, much 
current research continues to focus OjU one of the ear- 
liest studied candidate materials, Fe.cl Besides offering 
the possibility of room-temperature injection of electron 
spins, Fe has the potential advantage of forming a nearly 
lattice-matched epitaxial film on an important semicon- 
ductor substrate, GaAs. However, despite this early 
promise and after considerable research investment, mea- 
sured spin-injection efficiencies for Fe/GaAa_remain frus- 
tratingly low, typically no larger than 1%.lII3 The reason 
for this low efficiency is not definitively known, and much 
controversy surrounds its origin. A crucial issue, not yet 
settled, is whether the measured efficiencies reflect an in- 



trinsic upper limit, or are simply due to technical limita- 
tions that may be overcome or circumvented. For exam- 
ple, much early effort focused on the nature of interface 
layers inferred to be "magnetically dead", in the sense 
that they do not to CDJitribute to the total magnetic mo- 
ment of the Fe film.Elli3 The appearance of dead layers 
is consistent with the formation (thermodynamically fa- 
vorable in the bulk) of nonmagnetic FcAs complexes at 
an interface with the As-rich GaAs(001)-(2x4) surface. 
This reasoning led to the development of two strategies 
for suppressing As diffusion into the Fe film: passivatioa. 
of the As-rich surface by a surfactant such as sulfur;til 
and growing on the Ga-ric/i GaAs(001)-(4x6Vsiirface (or, 
equivalently. As decapping prior to growth) .113 Both ap- 
proaches lead to films with magnetization onset in the 
range 4-8 monolayers and with_.essentially the full mo- 
ment per Fe atom in all layers.EHnlj 

Notwithstanding the successful elimination of dead lay- 
ers, measured spin-injection efficiencies have, until re- 
cently, remained at or below the 1% level. Recently, 
Schmidt et al. have argued that a more fundamental 
limitation exists for spin injjeetion from a ferromagnetic 
metal into a semiconductor .Ej They have shown that in 
the purely diffusive regime (where spins are scattered 
much less frequently than electrons) the spin-injection 
efficiency from a ferromagnet (fm) into a semiconductor 
(sc) is proportional to the ratio of their conductivities, 
Csc/ctm- For ferromagnetic metals this ratio is of order 
lO^"*, and for typical device geometries suggests maxi- 
mum injection efficiencies of 1% or less. 

There remain several possibilities for circumventing 
this limitation on Fe sources. The first is to operate 
in the ballistic regime, where the contact resistance due 
to elastic backscattering at the interface (Sharvin re- 
sistance) will generally j«duce the metal-semiconductor 
conductivity mismatch.Ej Tang et al. have used concepts 
from mesoscopic transport to model injection into a two- 
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dimensional electron gas, and find clear evidence for bal- 
listic spinptsansport that otherwise vanishes in the diffu- 
sive limit 

A second possibility, relevant to the diffusive regime, 
is the use of tunnel contacts at the Fe/GaAs inter- 
face, which are expected to substantially reduce the con- 
ductivity i^ismatch and thereby increase the injection 
efhciency.U Related research avenues concern the role 
of intrinsic Schottky barriers in controlling the spin- 
dependent tunneling through an interface.Ej Indeed, Zhu 
et al. have recently demonstrated injection from Fe into 
GaAs with an efficiency of about 2% at room temper- 
ature; they suggest that the Schottky barrier formed 
between Fe and GaAs leads to a tunnel contact, thus, 
circumventing the problems of conductance mismatch.t2l 
Very recently, Hanbicki et al. have used a Schottky con- 
tact to an AlGaAs overlayer to inject spin into GaAs with 
efficiencies o£-30% at low temperature and 9% at room 
temperatureE3 

Since the prospects and limitations for spin injection 
from Fe into GaAs remain uncertain, we believe that fur- 
ther progress may benefit from a first-principles theoret- 
ical description of the interface. Our focus will be on the 
atomic structure of the interface and its resultant mag- 
netic character, especially during the first few monolayers 
of growth. Theoretical studies of buried interfaces are no- 
toriously difficult, for several reasons. First, experimen- 
tal probes can provide only indirect information about 
physical and electronic structure, and hence are of lim- 
ited utility for guiding theories. Moreover, real interfaces 
may — even after careful annealing — have atomic geome- 
tries very different from the ground state, hence the inter- 
face structure may depend on the precise growth history. 
In principle, one approach to this dilemma would be first- 
principles finite-temperature molecular dynamics simula- 
tion of film growth, which would properly account for the 
roles of deposition rate, surface diffusion, and incorpora- 
tion into the substrate. In practice, molecular dynamics 
using density functional theory can at best simulate pro- 
cesses for '--^lOO picoseconds — many orders of magnitude 
short of the experimentally relevant time scales, which 
may be milliseconds or longer. Kratzer and Scheffler have 
recently addressed this probleaji using a "first-principles 
kinetic Monte Carlo" methodO However, the complexity 
of applying this method to the growth of Fe on GaAs — 
with three atomic species and many different possible 
processes — makes such an approach not yet feasible. 

In this paper we approach the problem from two differ- 
ent limits. In Section || we consider the initial stages of 
interface formation, beginning with adsorption of isolated 
Fe adatoms on a bare GaAs substrate. We focus here on 
the potential-energy surface governing surface diffusion, 
and show that Fe-As chemistry may play a decisive role 
in the submonolayer regime. In Section III we consider a 
different limit, namely the thermodynamic ground-state 
interface structure of Fe films several monolayers thick. 
Here we concentrate on the magnetic character of the 
interface as a function of film thickness, and propose a 




FIG. 1. Potential energy surface for Fe on GaAs(OOl), 
calculated using GGA without spin polarization. The labeled 
points correspond to the configurations in Fig. Q. 



mechanism that accounts for the experimentally observed 
delayed ons et .. o£ f prromagnetic order for films of a few 
monolayers. tHloE3 In Section we discuss ways to rec- 
oncile our results from these two limiting regimes, and 
we propose future directions for further study. 



II. INITIAL STAGES OF GROWTH 

Our goal in this section is to identify and quantify those 
structural, magnetic, and chemical features that may be 
generic to the growth of Fe on GaAs. We do not attempt 
a definitive treatment of growth on a particular GaAs re- 
construction or at a particular temperature. Instead, we 
focus on the following four questions: (1) Does the initial 
adsorption of atomic Fe prefer metallic bonding at highly 
coordinated sites, or does covalent bonding to Ga or As 
prevail? (2) Does adsorbed Fe nucleate the formation of 
compact islands, or do adsorbate-substrate interactions 
dominate the growth? (3) Is surface diffusion of Fe likely 
to be significant at growth temperatures? (3) Are the 
Fe magnetic moments in the submonolayer regime gov- 
erned by strong Hund coupling (atomic moments) or by 
itinerancy effects (bulk moments)? 

We begin by studying Fe adsorption on a fictitious, but 
chemically reasonable, surface of GaAs(OOl): a (2 x 1) 
reconstruction consisting of bulk GaAs terminated by a 
dimerized Ga layer. This fictitious surface is analogous to 
the dimerized /3(2 x 4) and /32(2 x 4) As-rich surfaces, but 
is quite different from the more complex C(4 x 2) surface 
believed to be the basis for the c(|§-|X 2) reconstruction 
observed under Ga-rich conditions B3 

We use density-functional theory together with the 
generalized gradient approxiraation (GGA) for the 
exchange-correlation functional.^ For the calculations in 
this section we used TrouUier-Martins pseudopotentials 
and a plane-wave basis with a kinetic- energy cutoff of 
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FIG. 2. Relaxed configurations for Fe adsorbed at the four 
locations marked in Fig. ry. The gray and light gray spheres 
represent Ga and As atoms, and the smaller dark gray spheres 
represent Fe atoms. Highly strained bonds (in the range 
15-20% longer than bulk) are shown as dotted lines. The 
shaded area marks the region for which the potential-energy 
surface is plotted in Fig. hi. 

50 Ry, as implemented in the fhi96md codcil The sam- 
pUng of the surface Brillouin zone was equivalent to using 
64 fc-points for a (1 x 1) surface unit cell. 

To construct the potential-energy surface (PES) for ad- 
sorption of atomic Fe on the clean surface, we computed 
the total energy, Et, as a function of the adsorbate po- 
sition (X, Y) within the surface unit cell. To minimize 
interactions of the adsorbate with its periodic images we 
used a (2 X 2) supercell. For each adsorbate position 
{X,Y), we fully relaxed the Z coordinate of the adsor- 
bate and the positions of all Ga and As atoms in the 
top two layers of the surface. This procedure was re- 
peated for adsorbate positions sufficient to sample the 
PES with a resolution of about 0.25 A. For many adsor- 
bate positions, geometries with subsurface adsorption or 
with atomic positions exchanged were considered as well; 
in each case, the lowest-energy configuration was used to 
define Et{X,Y). 

In order to illuminate the detailed role — if any — of 
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FIG. 3. Relative total energy for Fe constrained to 
approach the Ga-dimer bridge site from directly above, as 
shown. Energies are with respect to the spin-unpolarized en- 
ergy for Fe far above the GaAs surface. The vertical distance 
between the two curves is the magnetic energy gain. With- 
out spin polarization there is a local barrier of ~0.4 eV for 
breaking the Ga dimer; when spin polarization is included 
this barrier is removed. 

magnetism in the adsorption energetics, we begin by 
first calculating the PES without allowing for spin- 
polarization. We then recalculate the PES while allowing 
for spin-polarization, and examine the differences that 
arise. If the adsorbate moment were completely local- 
ized and did not interact with the substrate, these two en- 
ergy surfaces would be identical except for an overall shift 
of the energies. Thus, the differences that we find (de- 
scribed below) directly refiect site-dependent magnetic 
interactions between the Fe moment and the GaAs sub- 
strate. 

Fig. |l| shows the PES calculated without spin polar- 
ization, for a portion of the surface unit cell shown in 
Fig. Ij. We restrict the plot to the vicinity of the dimer 
rows because the energy landscape between dimer rows 
is quite flat and considerably higher in energy. Along the 
dimer rows the PES is periodic with repeat length a/V2; 
Fig. shows slightly more than one full period along this 
direction. 

Each point on this PES corresponds to a different re- 
laxed geometry, determined solely by the in-plane adsor- 
bate position {X, Y). We focus on the four points marked 
in Fig. Il^, which correspond to the four geometries shown 
in Fig. 0. Point (a), in the center of the PES, corresponds 
to adsorption at the high-symmetry pedestal site, mid- 
way between two dimers along a dimer row. Although 
this site is highly coordinated, leading to four Fe-Ga and 
two Fe-As bonds of nearly equal length (2.5 A), it is ener- 
getically unfavorable. Indeed, this site is not even locally 
metastable: energy is gained by moving the adsorbate in 
any direction away from point (a) . 

The global minimum of the PES is a point (d) near 
the corner of the region shown in Fig. 0. This is 0.8 eV 
below point (a) and corresponds to the configuration of 
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FIG. 4. Potential energy surface for Fe on GaAs(OOl), 
calculated with spin polarization. The labeled points are the 
same as in Fig. nL and correspond very closely to the config- 
urations in Fig. 0. 

Fig. ||(d) , in which the Fe adsorbate has partially "kicked 
out" one Ga atom from a surface dimer and taken its 
place. The large energy gain from this process suggests 
that the formation of such Fe-Ga heterodimers may act 
as a strong local trap for Fe, strongly suppressing surface 
diffusion of the Fe adatoms. If the kicked-out Ga atom 
subsequently diffuses away to a more stable adsorption 
site, this trapping effect will be further enhanced. Such 
effects, if not kinetically hindered, will bias the growth 
toward nucleating many small Fe islands. 

The issue of kinetic barriers to forming Fe-Ga het- 
erodimers can be analyzed by considering the possible 
routes leading to point (d) along the PES. For example, 
the reaction pathway from point (a) to (d) proceeds via 
point (c), which corresponds to the configuration shown 
in Fig. |^(c). This is the transition state (the highest en- 
ergy configuration along the minimum-energy pathway) 
for the reaction taking (a) into (d), and is only 0.25 eV 
higher than the valley floor near (a). Part of this bar- 
rier arises from breaking the original Ga-Ga dimer bond, 
but apparently most of this energy cost is recovered by 
forming a more stable Fe-Ga heterodimer bond. Another 
1.0 eV is recovered at point (d) by forming, in addition, 
a new Fe-As backbond to the substrate. The kicked-out 
Ga atom may remain bonded to the Fe atom, but it is 
energetically much more favorable to break this Fe-Ga 
bond (at a cost of 0.6 eV) and diffuse to a more favor- 
able binding site — for a net gain up to 2.0 eV per Ga 
atom for the limiting case of incorporation into bulk Ga. 
Comparing these energy changes upon bond breaking, we 
conclude that Fe-As bonds are considerably more stable 
than Fe-Ga bonds. 

The final adsorbate site we discuss is point (b) , halfway 
between the pedestal site and the dimer bridge site. At 
this point, the Ga-Ga dimer bond has already broken to 
allow for more favorable Fe-Ga bonds to begin forming. 



FIG. 5. Comparison of the (a) spin-unpolarized and (b) 
spin-polarized valence electron density near the Fe atom for 
the configuration shown in Fig. ^(b). The contours are loga- 
rithmically spaced; lighter contours represent higher density, 
and are truncated near the Fe atom. Projected atomic posi- 
tions are marked by black circles; their size indicates proxim- 
ity to the plotting plane. 

Without spin polarization, the energy of this configura- 
tion is a local maximum, and represents a strong barrier 
to breaking Ga-Ga dimer bonds head-on. A similar bar- 
rier is found for a Fe adatom artificially constrained to 
approach the Ga-Ga dimer from directly above, as shown 
in in Fig. |3[ For this adsorption route from above, an 
energy cost of ~0.4 eV must be paid before the Ga-Ga 
dimer bond finally breaks; once paid, stable Fe adsorp- 
tion can occur at the bridge site. 

We turn now to the differences in the PES that arise 
from including spin polarization. We have repeated the 
calculation of the entire PES using spin-polarized GGA, 
allowing for any changes in the relaxed geometries. In 
general, we find only negligible changes to the relaxed 
atomic positions, and so the geometries shown in Fig. ^ 
continue to correspond to points on the new PES, shown 
in Fig. ^ Comparing the spin-unpolarized and spin- 
polarized energy landscapes (Figs. 1 and 4), several fea- 
tures deserve comment. The global minimum is the same 
for both, corresponding the Fe-Ga heterodimer formation 
via Ga "kick-out." With spin polarization, the pedestal 
site is again unstable against adsorbate motion in any 
direction. 

Differences are also apparent. The spin-polarized en- 
ergy surface appears more corrugated; this is due al- 
most entirely to a strong reduction of energies near point 
(b) (a local maximum without spin polarization) and 
the nearby dimer bridge site (a saddle point without 
spin polarization, but part of a low-energy trench with 
it). Thus by including spin-polarization, the barrier to 
head-on Ga-Ga dimer breaking is reduced to zero. The 
same difference is observed in dimer breaking by a Fe 
atom artificially constrained above the dimer: the barrier 
in Fig. ^ is eliminated by including spin polarization — 
demonstrating that this effect is not limited to surface 
diffusion. 

To understand why allowing for spin-polarization 
changes some parts of the PES but not others, we con- 
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FIG. 6. Relative total energies for a number of configura- 
tions used to define the potential-energy surfaces in Figs. 1 
and 4. The configurations within the bracket are all located 
near point (b) in Figs. 1 and 4 (see text for discussion). 

sider two plausible explanations: chemical effects and 
magnetic effects. By the former we mean contributions 
to the total energy that depend primarily on the total 
valence electron density; by the latter we mean contribu- 
tions related to the electron spin density. We focus on 
the configuration of Fig. |^(b), for which the change in 
the PES is particularly dramatic. In Fig. ^ we compare 
the spin-unpolarized and spin-polarized total valence- 
electron density in a plane containing the Fe adsorbate 
and its As and Ga neighbors. In both plots, Fe-As and 
Fe-Ga bonds are clearly visible, and the loss of both the 
Ga-Ga dimer bond and the Ga-As backbond is obvious. 
Most importantly, there are no large differences between 
the spin-polarized and unpolarized electron densities in 
this plane. We conclude that the changes in the PES due 
to spin-polarization cannot be attributed to changes in 
the valence electron density related to chemical bonds. 

To evaluate the role of magnetic effects in determin- 
ing the shape of the PES, we first consider how the total 
magnetic moment per cell varies with the Fe adsorbate 
position. In Fig. ^ we show the magnetic energy (the 
difference between spin-polarized and unpolarized ener- 
gies), as a function of total magnetic moment, for about 
50 different adsorption sites on the PES. Most of the 
sites have magnetic moments between 2 and 3 fiB , giving 
rise to magnetic energies between 0.4 and 0.8 eV. About 
ten sites have considerable higher moments, between 3.2 
and 3.4 /is, and correspondingly higher magnetic ener- 
gies, between 1.0 and 1.6 eV. Thus these sites, all located 
near point (b) on the PES, give rise to large differences 
in the spin-polarized and unpolarized energy landscapes. 
These differences can be attributed to the development 
of unusually large magnetic moments for adsorption sites 
in the vicinity of energy barriers in the spin-unpolarized 
PES; these sites generally involve highly strained or par- 
tially broken bonds. The same trend is observed for the 
constrained adsorption from above (Fig.^). The mag- 
netic energy gain at the configuration corresponding to 
the barrier is 2.5 eV, nearly twice as large as at the equi- 
librium adsorbate height. This difference again arises 




FIG. 7. (a) Potential energy surface for Fe on 

Fe/GaAs(001), calculated with spin polarization, (b) The 
initial preadsorbed surface structure, consisting of Fe-Ga het- 
erodimers alternating with Ga dimers. (c) Lowest-energy con- 
figuration for 1/2 monolayer Fe coverage, corresponding to 
the second Fe adsorbate at point (h) in the potential-energy 
surface. The shaded area marks the region plotted in (a). 

from the difference in magnetic moments: 3.9 fiB at the 
barrier configuration vs. 3.2 at equilibrium. Thus 
we conclude that magnetic effects play an important role 
in determining those parts of the PES corresponding to 
highly strained or partially broken bonds. 

We end this section by turning briefly to another pro- 
cess important during growth: the adsorption of Fe in 
the presence of preadsorbed Fe. We consider 1/4 mono- 
layer of Fe in its lowest-energy configuration, i.e. in- 
corporated as Fe-Ga heterodimers with (2x2) periodic- 
ity. To simplify the discussion, we assume a starting 
surface from which the kicked-out Ga atoms have de- 
tached and diffused away. We then recalculate the (spin- 
polarized) potential-energy surface for this preadsorbed 
surface, again fully accounting for geometrical relaxation. 

A portion of this PES is shown in Fig. |^(a), along with 
the initial preadsorbed surface structure in Fig. 0(b) . The 
area shown is roughly twice that of Fig. ^, and a general 
comparison of the corresponding portions reveals the fol- 
lowing. First, the overall corrugation of the preadsorbed 
PES is smaller than for the clean surface, i.e. preadsorbed 
Fe lowers the barriers to surface diffusion. The locations 
on this PES labeled (e) through (h) correspond to global 
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FIG. 8. Interface structures for the As-terminated 
Fe/GaAs(001) interface: Ideal (upper panel) and relaxed 
structures (lower panel). The gray and light gray spheres 
represent Ga and As atoms, and the dark gray spheres rep- 
resent Fe atoms. Highly strained bonds (15-20% longer than 
ideal) are shown as dotted lines. 

minima on the PES of the clean surface. Point (e), lo- 
cated at the position of the preadsorbed Fe, is here a lo- 
cal energy maximum. This suggests that Fe-Fe bonding 
is not yet favorable in this low coverage regime, at least 
relative to the further formation of Fe-Ga heterodimers. 
Points (f) through (h) are local minima, indicating that 
formation of Fe-Ga heterodimers remains favorable for 
the preadsorbed surface. The global minimum, point (h), 
corresponds to the staggered arrangement of Fe-Ga het- 
erodimers shown in Fig. ^c). Occupying all of these ad- 
sorption sites results in an Fe coverage of 1/2 monolayer. 

Although we do not explicitly calculate energy surfaces 
for higher coverages, we can make plausible inferences for 
their basic features based on the PES of Fig. |^. Even for 
a surface completely terminated by Fe-Ga heterodimers, 
the strong preference for Fe to form backbonds to As — as 
evidenced by the local minimum labeled (f) in Fig.|(a)— 
should persist. Thus we speculate that continued depo- 
sition of Fe will ultimately lead to Fe-Fc dimers atop 
the As-terminated substrate. An obvious consequence of 
such a configuration would be the release of Ga (a full 
monolayer in this surface model). Moreover, in Section 
III we will show that for Fe films of several monolayers, 
a floating Ga (or As) adlayer can act as a surfactant, 
lowering the surface energy by as much as 2 eV/(lxl) 
cell. 



To conclude this section, we have shown that the initial 
stages of adsorption of Fe are dominated by strong local 
Fe-As chemistry. For a surface terminated by Ga, this 
chemistry leads to facile breaking of surface Ga dimers 
by sequential kicking-out of the two Ga atoms, first form- 
ing Fe-Ga heterodimers and, finally, surface Fe dimers 
bonded to subsurface As. Excess Ga may be released 
as a result of this rebonding mechanism. Magnetism 
plays an auxiliary role in the process by lowering the 
potential-energy barriers to breaking apart surface Ga 
dimers. Since the formation of Fe-As bonds leads to effi- 
cient trapping of Fe — especially if the excess Ga diffuses 
away — surface diffusion of Fc will probably be strongly 
suppressed. 

III. GROUND-STATE INTERFACE STRUCTURE 

As a complementary approach to the study of single Fe 
adsorption and diffusion, we have studied epitaxial inter- 
faces of Fe/GaAs. In this section, our focus is threefold: 
(1) to determine the stability and magnetic character of 
Fe films for different interface structures and film thick- 
ness; (2) to investigate the role of As or Ga adlayers on 
Fe films, in particular the extent to which they may ac- 
count for the observed out-diffusion of substrate atoms 
to the surface; (3) to examine the possibility of antifer- 
romagnetic order as the origin of the observed magnetic 
quenching of Fe films with thicknesses of just a few mono- 
layers. 

The Fe/GaAs interface structures considered in this 
section are shown in Fig. ^ for the the As-terminated 
GaAs(OOl) interface. Model A is an atomically abrupt in- 
terface of bcc Fe and zinc-blende GaAs. Because the lat- 
tice constant of bcc Fe (a = 2.866 A) is almost half that 
of the substrate (a = 5.654 A), the epitaxial relationship 
is Fe(001)(100) || GaAs(001)(100) and the atomic density 
of each Fe layer is twice that of the substrate layer, with 
a strain of only 1.3%. With respect to the bcc Fe lat- 
tice, there are vacancy sites in the adjacent GaAs lattice. 
Models B and C are built by filling these vacancy sites 
with Fe atoms one by one; we will refer to these inter- 
faces as "partially intermixed" and "fully intermixed," 
respectively. For each of these models, we considered 
film thicknesses ranging from 0.5 ML [one Fe atom per 
(1 X 1) of GaAs(OOl)] up to 3.5 ML. We also considered 
three analogous models for the Ga-terminated interface. 
Both As- and Ga-terminated interface models may have 
relevance for experiments with the more commonly used 
Ga-rich surfaces: for example, we have already seen in the 
previous section that substitutional displacement of Ga 
on a Ga-rich surface leads to an As-terminated Fe/GaAs 
interface. 

For the calculations in this section we again used the 
generalized gradient approximation (GGA), here with an 
ultrasoft pseudopotentialEEl for Fe. The electronic wave 
functions and densities were described by a planewave 



6 




.9 sH Ga-terminated 

CO 

E 
o 



5 4- 




2- 



12 3 4 

Fe coverage (ML) 

FIG. 9. Formation energies [in eV/(lxl)-cell] of Fe films 
on tlie GaAs(OOI) surface given as a function of interface 
structure (A,B,C) and Fe coverage. The formation energies 
given are taken at the center of the thermodynamically al- 
lowed range of the As chemical potential: in the As-rich limit 
condition the formation energies decrease (increase) by 0.16 
eV for the As-terminated (Ga-terminated) Fe/GaAs interface. 



basis with cutoff energies of 16 Ry and 160 Ry, respec- 
tively. For each structural model, six or seven atomic 
layers were used for the substrate; the bottom As or Ga 
layer was passivated by pseudohydrogen atoms and fixed 
during the structural relaxation. The relative stability of 
each slab is given by its formation energy, 
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where Et is the total energy of the slab, Ni the number 
of atoms of each chemical type, and fii their chemical 
potentials. To eliminate the contribution of the pseu- 
dohydrogen layer to the formation energy, we define the 
formation energy of each structural model with respect 
to a common reference structure. 
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The reference structure is an ideal GaAs slab passivated 
on both sides (hence the factor 1/2) by pseudohydrogen 
atoms. Assuming the thermodynamic equilibrium con- 
dition, UGaAs = ^J■As + MGa, and taking /XCaAs and /Xpe 

from bulk structures, the formation energy is expressed 
as a function of fi^s only within the thermodynamically 
allowed range ^caAs - MGa(buik) < A^As < MAs(buik)- Here, 
the lower (upper) limit corresponds to the Ga-rich (As- 
rich) environment, and /iGa(buik) and /iAs(buik) are again 
determined from their bulk structures. 

Figure || shows the calculated formation energies as a 
function of Fe film thickness. We first discuss Fe films 
with the As-terminated interface. At a Fe coverage of 
0.5 ML model C is most stable, while models A and B 
are ~1 eV higher. This energetic ordering persists up to 
1 ML coverage, but changes above that point: the forma- 
tion energy of model C increases by 1.8 eV from 1 ML to 
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FIG. 10. Formation energies and magnetic moments for 
four selected adlayer structures. The Fe film thickness, sub- 
strate termination, and interface model is listed underneath 
each structure. Formation energies and magnetic moments 
are shown for bare Fe films and for As- and Ga-adlayers. The 
formation energies are calculated at the center of the ther- 
modynamically allowed range of the As chemical potential 
{fj. As)- The changes at the limiting values of ^as are indi- 
cated by small circles (filled circles correspond to the Ga-rich 
limit, open circles to the As-rich limit). 



2 ML, while that of the A and B become nearly indepen- 
dent of coverage. After the formation of 2 ML, the for- 
mation energy of all three models does not change much 
by adding more Fe layers. Thus, the low-energy interface 
structure turns out to be model A — the abrupt interface. 
For Fe films with the Ga-terminated interface, a similar 
trend apphes but with a different outcome. At 0.5 ML, 
models B and C are equally stable, and at 1 ML model 
C is in fact the most stable. With additional Fe layers, 
model C becomes very unstable, similar to the case of the 
As-terminated interface. At higher coverages, however, 
model B — the partially mixed interface — is most stable. 

These results can be understood as a competition be- 
tween maximizing the coordination of Fe atoms and min- 
imizing the concentration of atoms in the interfacial re- 
gion. At low Fe coverage (< 1 ML), model C is energeti- 
cally most favored, because Fe atoms can maximize their 
coordination with substrate atoms. At higher coverages, 
where the interface is well defined, this arrangement be- 
comes unstable relative to the less intermixed interfaces. 
The reason for this crossover is the two extra Fe atoms 
per (1x1) surface cell (relative to model A). The ex- 
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cess electrons from the extra Fe atoms fill antibonding 
orbitals and thus weaken the interface bonding. This 
weakened bonding is also evident in the relaxed interface 
separations: for both As- and Ga-terminated interfaces, 
the position of the first Fc monolayer is ^0.8 A higher 
in model C and ^0.4 A higher in model B, compared 
to model A (see the relaxed structm^es in Fig. ||). Simi- 
lar reasoning explains why model B is more stable than 
model A for the Ga-terminated interface, while the op- 
posite holds for the As-terminated interface. 

The present results emphasize that the low-energy 
atomic structure at low Fe coverage may not be extrap- 
olated to high coverages. This should be true also in 
real growth situations, since Fe atoms at low coverages 
would take positions one or two layers deep so as to 
maintain maximal coordination. Further deposition of 
Fe leads to an partially mixed interface that ultimately 
becomes unstable. Therefore, it is likely that substan- 
tial rearrangement of the atomic structure occurs dur- 
ing the film growth, provided the temperature is suffi- 
ciently high. This may partly account for the experimen- 
tal observationnBcZlna that substrate atoms, especially As 
atoms, diffuse out to the surface during the Fe growth 
even for room-temperature deposition. 

To follow this reasoning, we examined the effect of a 
Ga or As adlayer on top of the Fe film. In Fig. |o| we 
compare the formation energies of four selected interface 
structures with and without an adlayer. Both Ga and 
As adlayers stabilize the surface substantially, regard- 
less of the interface structure and film thickness. This is 
not surprising, since adlayers increase the coordination 
of surface Fe atoms. The energy gain from As adlayers is 
particularly large, amounting to more than 1.5 eV/(l x 1) 
in most cases. Adlayers on 1 ML of Fe show strong de- 
pendence on the adatom site: the low-energy site is the 
position extended from the substrate GaAs lattice, sug- 
gesting a strong covalent bonding between the substrate 
and the adlayer, mediated by the intervening Fe layer. 

The results for adlayer structures, together with the 
interface energetics, provide strong theoretical evidence 
that substitutional adsorption and/or atomic exchange 
is an essential process during the Fe raowth — a finding 
which was anticipated by experimentsEj and which we 
explicitly demonstrated in Section By substitutional 
adsorption processes, Fe atoms can maximize their coor- 
dination at every stage of the growth. These processes re- 
sult in the segregation of substrate atoms to the surface, 
and simultaneously facilitate optimal interfacial atomic 
densities and thus relatively stable interface structures. 
Our calculations have shown that such processes are in- 
evitable and that the energy gain from them is quite sub- 
stantial. 

We now turn to magnetic properties of the interface 
structures. Calculated spin moments (per Fe atom) of 
various model structures are given in Figs. |l| and 
For bare Fe films without Ga or As adlayers, average Fe 
spin moments are much enhanced compared to the calcu- 
lated bulk value for Fe of 2.33 hb- Fig. |ll| shows that for 




12 3 4 

Fe coverage (ML) 

FIG. 11. Magnetic moments of Fe films on GaAs(OOl) as 
a function of interface type and Fe coverage. 

TABLE L Atom-resolved spin moments (in /ib) for 
Fe/GaAs interface structures. The layer number is given with 
respect to the top substrate layer. The induced spin moments 
of substrate atoms are given with a species label. 

Layer As:A 3.0 ML Ga:B 3.5 ML As:A 2.0 ML \ As 



3 


3.03, 3.00 


3.02, 2.99 


-0.08 (As 


2 


2.51, 2.38 


2.49, 2.43 


2.52, 2.42 


1 


2.54, 2.49 


2.49, 2.42 


2.40, 2.29 





-0.04 (As) 


2.72, -0.06 (Ga) 


-0.04 (As 


-1 


0.02 (Ga) 


0.00 (As) 


0.02 (Ga) 


-2 


0.01 (As) 


0.03 (Ga) 


0.02 (As) 


-3 


0.00 (Ga) 


0.01 (As) 


0.00 (Ga) 



films of 3 ML this enhancement is still substantial, and 
relatively insensitive to the specific interface structure or 
substrate termination. For Fe films with adlayers. Fig. |l^ 
shows that both As- and Ga-adlayers suppress total mag- 
netic moments, by as much as 1 /is in several cases. 

In order to examine local variations in the total mag- 
netic moments, we calculated partial moments within a 
sphere centered on each atom, in analogy to the muffin- 
tin sphere in all-electron approaches. We used a sphere 
radius of 1.2 A, which is slightly shorter than the min- 
imum bond length between atoms. Table | shows the 
results for three selected structures. Local Fe moments 
are significantly enhanced (~3.0 /is) at the surface layer. 
For the relatively thin Fe films considered in this work, 
the buried and interfacial Fe layers also have sizable en- 
hancement. The adlayer suppresses spin moments not 
only of the top-layer Fe atoms but also slightly of the 
second-layer Fe atoms. Both the enhancement and sup- 
pression can be understood in terms of changes in the Fe 
d-band width due to changes in coordination and sym- 
metry. On the other hand, substrate atoms have small 
induced moments up to 3-4 layers deep, always with neg- 
ative spin moments for the interfacial atoms. The sign of 
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FIG. 12. Formation energies of ferromagnetic (FM) and 
antiferromagnetic (AF) Fe films in two different Fe/GaAs 
structures (As- and Ga-terminated abrupt interfaces). Three 
types of adlayers are considered (none, As, Ga). 



bimolecular units with two Fe atoms in the same plane, 
plus (me As and one Fe atom on each side of the Fe 
plancEJ Thus it is plausible that AF order in Fe/GaAs 
interfaces at low coverage arises from the formation of 
the Fe2As-like structures. This local similarity is broken 
by additional Fe adsorption, so that ferromagnetic order 
is ultimately favored. 

To conclude this section, we have shown that the low- 
energy interface structure is different for low and high 
Fe coverages. The addition of Ga- or, especially, As- 
adlayers substantially lowers the surface formation, pro- 
viding a theoretical basis for the experimentally observed 
out-diffusion of substrate atoms to the surface. We have 
shown that Fe films on the GaAs substrate usually as- 
sume a ferromagnetic ground state even at very low cov- 
erages, while for certain cases an As adlayer can induce 
antiferromagnetic order in the Fe film. Since the lat- 
ter cases generally have lower formation energies and are 
thus likely to form, we propose that the observed mag- 
netic quenching for very low Fe coverage may be due to 
the formation of antiferromagnetic order. 



spin moments changes to positive at deeper layers, sug- 
gesting the formation of a spin-density wave. 

Recent magnetic measurements have observed a de- 
layed onset of the magnetic phase at ^2 ML for Fe films 
grown 0|ft,the Ga-rich GaAs(001)-(2 x 6) and (4 x 2) 
surfaces.Ea Of the various surface structures we have con- 
sidered here, one shows a similar quenching effect. The 
As-terminated interface with 1 ML of Fe (the first struc- 
ture in Fig. p^ ) has a zero net moment when capped by an 
As adlayer. Since this structure also has the lowest for- 
mation energy for an Fe coverage of 1 ML, it might be an- 
ticipated that the delayed onset is due to the quenching of 
magnetism from the strong covalent bonding between the 
As atoms via Fe atoms. However, there is another pos- 
sibility, namely the formation of antiferromagnetic (AF) 
order. We have considered the possibility of AF order 
for the simplest structures with 1 ML of Fe in the same 
plane, so that one Fe atom has spin up and the other 
down. The results, summarized in Fig. indicate that 
films with an As adlayer become more stable by 0.05- 
0.09 eV/(l X 1) upon the formation of AF order, while 
bare films and films with a Ga adlayer become unstable 
by 0.2-0.5 eV/(l x 1). We note that the lowest energy 
structure with quenched magnetism is unstable toward 
AF order. We propose on this basis that the observed 
delayed magnetic onset is due to the initial formation of 
AF order, and tightly correlated with the out-diffusion 
of As atoms to the surface. 

While the microscopic origins driving the formation of 
AF order are not clear in detail, it appears that the sur- 
face As atoms mediate the AF order between Fe atoms. 
A related observation is that the AF-stabilized surface 
structure is locally similar to the tetragonal Fe2As struc- 
ture in the AF ground state: the unit cell consists of 



IV. SUMMARY AND FUTURE DIRECTIONS 

We have taken two complementary approaches to un- 
derstanding the nucleation and growth of Fe on GaAs — 
one focusing on the behavior of single Fe adatoms de- 
posited on clean and partially preadsorbed GaAs(OOl), 
the other focusing on the interface structure of complete 
Fe films at coverages up to several monolayers. Although 
a detailed growth history of Fe/GaAs interfaces cannot 
yet be described, our studies suggest four generic features 
that may play an important role. 

First, for low coverages we have identified a very strong 
driving force for Fe to be highly coordinated. For sin- 
gle adatoms, this tendency is strong enough to sponta- 
neously break surface Ga-Ga and Ga-As bonds in order 
to form Fe-As bonds. Although magnetic effects do not 
dominate this chemistry, they play an interesting auxil- 
iary role by essentially removing energetic barriers from 
the reaction pathways. For half-monolayer films con- 
strained by symmetry, these local surface chemical re- 
actions are not possible, and the Fe atoms respond by 
occupying subsurface sites with high coordination to ei- 
ther Ga or As. 

A second generic feature is the crossover from a pref- 
erence for strongly intermixed films to less intermixed or 
even abrupt films. This crossover occurs between one 
and two ML for both As- and Ga-terminated interfaces. 
Its origin is not magnetic: the films exhibit moments per 
atom larger than the bulk value even for very low cov- 
erages, and the moments converge to their bulk value 
long after the crossover occurs. Rather, it arises from 
the competition between maximizing the coordination of 
Fe atoms (which favors intermixing) and minimizing the 
amount of excess interfacial Fe (which favors abrupt in- 
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terfaces). For well-defined Fe films of 2 ML or more, the 
latter effect dominates and sharper interfaces become en- 
ergetically preferred. 

A third finding, common to all interfaces we have stud- 
ied, is that Ga- and As-adlayers dramatically reduce the 
formation energies of Fe films. This stabilizing effect 
occurs for both Ga- and As-terminated interfaces, for 
both intermixed and abrupt interfaces, and for all film 
thicknesses considered. It is especially striking for As- 
adlayers, which can reduce the film formation by as much 
as 50%. We also find that this stabilization is generally 
accompanied by a suppression of the total magnetic mo- 
ment of the film; since this is due to reduced local mo- 
ments in the topmost layer or two, the effect is largest 
for thinner films. Indeed, for 1-ML Fe films with an As- 
adlayer, antiferromagnetic order can be more stable than 
ferromagnetic order. 

Taken together, these three generic features collec- 
tively imply a fourth: the diffusion of Ga or As atoms 
from the interface to the surface of the Fe film. We have 
shown explicitly how Ga can be released by adsorption 
of individual Fe adatoms and the subsequent "kick-out" 
of Ga from surface dimers. Even for GaAs surfa|Ce re- 
constructions that do not consist of Ga dimers,E3 we 
speculate that the same strong Fe-As chemistry would 
again lead to the release of surface Ga atoms. We have 
also shown that during the growth — between one and two 
ML — a spontaneous rearrangement of the interface mor- 
phology is likely to occur, again leading to the release of 
either Ga or As. Although we do not speculate about the 
details of this atomic rearrangement, we have also shown 
that ultimately the liberated Ga or As is likely to play 
the role of a floating surfactant layer. 

Finally, we mention a possible avenue for further re- 
search. One difficulty with theoretical studies of inter- 
face is the paucity of macroscopic observables that can 
be directly related to the microscopic interface struc- 
ture. Schottky barriers are extremely sensitive probes 
of interface structure, varying-by as much 25% for local 
changes in interface geometry.tll Schottky barrier heights 
have been measured in Fe/GaAs interfaces to be of order 
0.7 eVj23 and thus represent a useful probe of interface 
microstructure. Moreover, Schottky barriers may actu- 
ally be a necessary ingredient for circumventing the in- 
trinsic limitations on spin injection from a ferromagnetic 
metal into a semiconductor. A theoretical understand- 
ing of their dependence on Fe/GaAs interface structure — 
including substrate termination and reconstruction, de- 
gree of intermixing, and magnetic character — would be a 
great asset. 

V. ACKNOWLEDGMENTS 

This work was funded in part by ONR and in part by 
the Deutsche Forschungsgemeinschaft. Computational 
work was supported in part by a grant of HPC time 



from the DoD Major Shared Resource Center ASCWP. 
S.C.E. and S.H.L. acknowledge generous support from 
the Alexander von Humboldt Foundation. 



* Permanent address: Computational Science and Engineer- 
ing Center, Samsung Advanced Institute of Technology, 
P.O. Box 111, Suwon 440-600, South Korea 

^ G. A. Prinz, Science 282, 1660 (1998) 

^ R. Fiederling, M. Keim, G. Reuscher, W. Ossau, G. 
Schmidt, A. Waag, and L. W. Molenkamp, Nature (Lon- 
don) 402, 787 (1999). 

^ B. T. Jonker, Y. D. Park, B. R. Bennett, H. D. Cheong, G. 
Kioseoglou, and A. Petrou, Phys. Rev. B 62, 8180 (2000). 

*Y. Ohno, D. K. Young, B. Beschoten, F. Matsukura, H. 
Ohno, and D. D. Awschalom, Nature (London) 402, 790 
(1999). 

5 C. Haas, IBM J. Res. Dev. 14, 282 (1970). 
® G. A. Prinz and J. J. Krebs, Appl. Phys. Lett. 39, 397 
(1981). 

^ P. R. Hammar, B. R. Bennett, M. J. Yang, and M. Johnson, 

Phys. Rev. Lett 83, 203 (1999). 
* Y. B. Xu, D. J. Freeland, E. T. M. Kernohan, W. Y. Lee, 

M. Tselepi, C. M. Guertler, C. A. F. Vaz, J. A. C. Bland, S. 

N. Holmes, N. K. Patel, and D. A. Ritchie, J. Appl. Phys. 

85, 5369 (1999). 
^ J. J. Krebs, B. T. Jonker, and G. A. Prinz, J. Appl. Phys. 

61, 2596 (1987). 

A. Filipe, A. Schuhl, and P. Galtier, Appl. Phys. Lett. 70, 
129 (1997). 

" G. W. Anderson, M. C. Hanf, and P. R. Norton, Phys. Rev. 
Lett 74, 2764 (1995). 

M. Zolfl, M. Brockmann, M. Kohler, S. Kreuzer, 
T. Schweinbock, S. Miethaner, F. Bensch, and G. 
Bayreutherp, J. Mag. Magn. Mat. 175, 16 (1997). 
" E. M. Kneedler and B. T. Jonker, J. Appl. Phys. 81, 4463 
(1997). 

" Y. B. Xu, E. T. M. Kernohan, D. J. Freeland, A. Ercole, M. 
Tselepi, and J. A. C. Bland, Phys. Rev. B 58, 890 (1998). 

G. Schmidt, D. Ferrand, L. W. Molenkamp, A.T. Filip, and 

B. J. van Wees, Phys. Rev. B 62, R4790 (2000). 

H. X. Tang, F. G. Monzon, R. Lifshitz, M. C. Cross, and 
M. L. Roukes, Phys. Rev. B 61, 4437 (2000). 

^''E.I. Rashba, Phys. Rev. B 62, R16267 (2000). 

A. Hirohata, Y. B. Xu, C. M. Guertler, and J. A. C. Bland, 
J. Appl. Phys. 87, 4670 (2000). 

H. J. Zhu, M. Ramsteiner, H. Kostial, M. Wassermeier, H.- 
P. Schonherr, and K. H. Ploog, Phys. Rev. Lett 87, 016601 
(2001). 

A. T. Hanbicki, B. T. Jonker, G. Itskos, G. Kioseoglou, and 

A. Petro, Appl. Phys. Lett. 80, 1240 (2002). 

P. Kratzer and M. Scheffler, Computing in Science and 
Engineering 3, 16 (2001). 

E. M. Kneedler, B. T. Jonker, P. M. Thibado, R.J. Wagner, 

B. V. Shanabrook, and L. J. Whitman, Phys. Rev. B 56, 
8163 (1997). 



10 



S. H. Lee, W. Moritz, and M. Scheffler, Phys. Rev. Lett 

85, 3890 (2000). 

J. P. Perdow, J. A. Chcvary, S. H. Vosko, K. A. Jackson, 
M. R. Pederson, D. J. Singh, and C. Fiolhais, Phys. Rev. 
B 46, 6671 (1992). 

M. Bockstedte, A. Kley, and M. Scheffler, Computer 
Physics Commun. 107, 187 (1997). 

D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

S. A. Chambers, F. Xu, H. W. Chen, L M. Vitomirov, S. 
B. Anderson, and J. H. Weaver, Phys. Rev. B 34, 6605 
(1986). 

T. L. Monchesky, B. Heinrich, R. Urban, K. Myrtle, M. 
Klaua, and J. Kirschner, Phys. Rev. B 60, 10 242 (1999). 
F. Bensch, G. Garreau, R. Moosbuhler, G. Bayreuther, and 

E. Beaurepairc, J. Appl. Phys. 89, 7133 (2001). 

R. W. G. Wyckoff, Crystal Structures, 2nd Ed., VoL 1, 
(Interscience Publishers, New York, 1963), p. 360. 
D. R. Heshnga, H. H. Weitering, D. P. van der Werf, T. M. 
Klapwijk, and T. Hibma, Phys. Rev. Lett 64, 1589 (1990). 
J. A. C. Bland, A. Hirohata, C. M. Guertler, Y. B. Xu, 
and M. Tselepi, J. Appl. Phys. 89, 6740 (2001). 



11 



